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Abstract: We study soft-gluon radiation in top quark decay within the framework 
of perturbative fragmentation functions. We present results for the 6-quark energy 
distribution, accounting for soft-gluon resummation to next-to-leading logarithmic 
accuracy in both the MS coefficient function and in the initial condition of the per- 
turbative fragmentation function. The results show a remarkable improvement and 
the 6-quark energy spectrum in top quark decay exhibits very little dependence on 
factorization and renormalization scales. We present some hadron-level results in 
both xb and moment space by including non-perturbative information determined 
from e + e~ data. 
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1. Introduction 

Heavy-flavour and in particular top quark physics is presently one of the main fields 
of investigation in theoretical and experimental particle physics. The current ex- 
periments at the Tevatron accelerator and, ultimately, at LHC [1] and e + e~ Linear 
Colliders [2] will produce large amounts of top quark pairs, which will allow one 
to perform improved measurements of top properties, such as its mass. For this 
purpose, precise calculations for top production and decay processes are mandatory. 

While fixed-order calculations reliably predict total cross sections or widths, 
differential distributions typically contain large logarithms associated with soft or 
collinear parton radiation. For heavy-quark production processes, although the quark 
mass m (much larger than the QCD scale A) acts as a regulator for the collinear sin- 
gularity, event shapes still contain large \n p (Q 2 /m 2 ) (with p < n) terms, Q being 
a typical scale of the process, which make fixed-order predictions unreliable when 
Q 3> m. Such logarithms can be resummed by using the perturbative fragmentation 
approach [3], which factorizes the rate of heavy-quark production into the convolution 
of a coefficient function, describing the emission of a massless parton, and a pertur- 
bative fragmentation function D(/iF,m), where \if is the factorization scale. The 
perturbative fragmentation function expresses the transition of the massless parton 
into the massive quark, and its value at any scale hf can be obtained by solving the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations [4,5] once 
an initial condition at a scale /i F is given. In ref. [3] the large \n(Q 2 /m 2 ) collinear 
logarithms were resummed to next-to-leading logarithmic (NLL) accuracy in the case 
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of heavy quark production in e + e~ collisions and an explicit next-to-leading order 
(NLO) expression for .D(//p,m), which was argued to be process independent, was 
given. More recently, ref. [6] has established the process independence in a more 
general way. 

The approach of perturbative fragmentation has been extensively used for e + e~ 
annihilation [6-10], hadron collisions [11,12], photoproduction [8,13] and, more re- 
cently, for bottom quark production in top quark decay t — > bW [14]. 

The initial condition of the perturbative fragmentation function and the coeffi- 
cient function, though free of the collinear large ln(Q 2 /m 2 ), contain large logarithms 
which are due to soft-gluon radiation. The ones contained in the initial condition 
are process independent [6], and were already included in the top-to-bottom decay 
process in [14], where the complete 0(a s ) calculation of the t — > bW(g) process was 
also performed. The large logarithms contained in the coefficient function are instead 
process dependent, and have to be evaluated for every specific case. It is precisely 
the purpose of this paper to extend the analysis of ref. [14] and to present results 
for soft-gluon resummation in the coefficient function. This will allow, together with 
the results of [6], to complete the evaluation of soft-gluon effects to NLL accuracy 
for the top-to-bottom decay process, and to investigate the impact on the 6-quark 
energy distribution. 

After the fragmentation of heavy quarks in e + e~ collisions considered in [6], this is 
the first process whose large logarithms (both collinear and soft) are fully resummed 
to NLL accuracy within the perturbative fragmentation function formalism. The 
consistency of this perturbative description with the one used in the e + e~ process 
makes it possible to fit non-perturbative information from e + e~ data and use it to 
make predictions. We shall therefore be able to predict the spectrum for 6-flavoured 
hadron energy distributions in top decay using e + e~ experimental data from LEP. 

The outline of the paper is as follows. In section 2 we review bottom quark 
production in top quark decay within the framework of perturbative fragmentation. 
In section 3 we present analytic results for the NLL soft-gluon resummation of the 
coefficient function in top decay. At the end of the section we also comment on 
the relation between our results and previous work on soft-gluon resummation in 
heavy-flavour decay [15-18]. In section 4 we show the 6-quark energy distribution 
in top decay and investigate the impact of soft-gluon resummation. In section 5 
we discuss inclusion of non-perturbative effects and present results for 6-flavoured 
hadron spectra in top decay. In section 6 we summarize our main results. 

2. Perturbative fragmentation and top quark decay 

We consider top decay into a bottom quark and a real W boson plus, to order as, a 
gluon: 

t(pt)^b(p b )W(p w )(g( Pg )) (2.1) 
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and define the bottom and gluon normalized energy fractions x b and x 
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^ f 1 dz \ 1 dTi \1 n Xb 

{x b ,m t ,m w ,m b ) = > / — - — — {z,m t ,m w ,fj, F ) Di[—,n F ,m b 
j Jx b z T B dz \ z 



1 — w mf 1 — w mf 

where u> = m\ v jm\. Neglecting the b mass, we have < x b>g < 1. 

Order as corrections to the decay process (2.1) were considered in [14]. It was 
observed there that, since m b <C m t , one can readily neglect m b /m t power suppressed 
terms, but on the other hand it is important to resum to all orders terms enhanced, 
at order as, by the presence of ln(m 2 /m b ). Such a resummation was performed 
in [14] by employing the perturbative fragmentation formalism [3,6]: The differential 
width for the production of a massive b quark in top decay is written in terms of the 
the convolution 

1 dT , . ^ f 1 dz T 1 dTi 

( 

T B dx b 

+0((m b /m t Y) , (2.3) 

where T B is the width of the Born process t — > bW, dfi/dz is the differential 
width for the production of a massless parton i in top decay with energy fraction z, 
Di(x, fj,F,iTib) is the perturbative fragmentation function for a parton i to fragment 
into a massive b quark, fip is the factorization scale. The term O ((m b /mt) p ) on the 
right-hand side stands for contributions that are suppressed by some power p (p > 1) 
of m b in the m b m t regime. Of course, non-perturbative corrections of the type 
A/m t and A/m b are understood on the right-hand side of eq. (2.3). We shall use 
everywhere a branching fraction B(t — > bW) = 1, and only include i = b in the above 
summation. 

The massless differential distribution (\/Y B ) dTi/dz (which is what we shall also 
refer to as "coefficient function") is defined in the MS factorization scheme after 
subtraction of the collinear singularities. It has been calculated at order as in [14]. 
In the following we shall often use its Mellin moments, defined by 

T N = f dz z N ~i-Lf(z). (2.4) 
Jo 1 B dz 

In moment space the convolution (2.3) can then be rewritten as 

rjv(m t , m w , m b ) = t N (m t , m w , /j, F )D biN (fi F , m b ) . (2.5) 

The perturbative fragmentation function D b (x, /if, mb) at any scale fi F can be 
obtained by solving the DGLAP equations. As shown in [6], as long as one can 
neglect contributions proportional to powers of (m b /m t ) p , the initial condition for 
the perturbative fragmentation function, which we evaluate at a scale //of, is process 
independent (but scheme dependent). In the MS scheme it reads [3]: 



D b ni (x,HoF,m b ) = 5(l-x)^ 



2tt 



1 + x 2 ( fJ* F 



l n ^_21n(l-x)-l 



x \ mt 



■ (2-6) 

+ 
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The solution of the DGLAP equations in the non-singlet sector, for the evolution 
from the scale /Iqf to /If, is given in moment space by: 



D btN (fj, F , m b ) = D^hof, m b ) exp 



3 (o) 

N 



In 



+ 



Q=s(/4)f) 

27T&0 " tts(/4) 
p (l) 27T&! (0) 



(2.7) 



In eq. (2.7) and are the Mellin transforms of the leading and next-to-leading 
order Altarelli-Parisi splitting vertices, and their explicit expression can be found, 
e.g., in [3]. b and bi are the first two coefficients of the QCD /5-function 



33-2n f , 153 -19n* 



(2.J 



12tt ' _± 24tt 2 ' 

which enter the following expression for the strong coupling constant at a scale Q 2 : 

1 [ 6i In [ln(Q 2 /A 2 )] 1 



& hi(Q 2 /A 2 



bl ln(Q 2 /A 2 



(2.9) 



Equation (2.7) resums to all order terms containing large \n(/i F / iXq F ) ■ I* 1 particu- 
lar, leading (agin" (fxp / ' Hq F )) and next-to-leading (ag In™ -1 (fx F / lxq F )) logarithms are 
resummed. Setting, as done in [14], jip — tth an d I^of — "fy>, one resums the large 
In (m 2 /m 2 ) terms with NLL accuracy. These are indeed the large collinear logarithms 
exhibited by the fixed-order calculation with a massive b quark [14]. 



3. Soft-gluon resummation 

In this section we address the problem of soft-gluon resummation in top quark de- 
cay. The MS coefficient function computed in [14] and the initial condition of the 
perturbative fragmentation function (2.6) contain terms behaving like 1/(1 — x)+ or 
[ln(l — x)/(l — x)] + , which become arbitrarily large when x approaches one. This is 
equivalent to contributions proportional to In N and In 2 N in moment space, as can 
be seen by writing the MS coefficient function [14] in the large- N limit 1 : 



r N (m t ,m w ,n F ) 



^{ 2 ln 2 iV + 
2tt I 



4 7e + 2 - 4 ln(l - w) - 2 In 

Hf 



+ K(m t) m w , ix F ) + O 




XnN 
(3.1) 



following [18], we note that, by denning n = Nexp^s), we could rewrite this expression in 
terms of ln(n) rather than In TV, with no terms explicitly appearing. 
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where 7^ = 0.577... is the Euler constant and w = m^/m 2 , as defined in section 1. 
In eq. (3.1) we have introduced the function K(m t ,mw, (J>f), which contains terms 
which are constant with respect to N. It reads: 

,2 



K(m t ,mw,fJ>F) 



3 \ TYl ^ 

- - 2 7jr ln-£ + 2 7 | + 2 lE [1 - 21n(l - w)\ 
z J Hp 



+ 21nwln(l -w)-2 
+ 4Li 2 (l 



w 



l + 2w 

-2 



ln(l — w) — 



2w 
1 — w 



Inw 



w 



(3.2) 



The x — > 1 (AT — > 00) limit corresponds to soft-gluon radiation in top decay. 
These soft logarithms need to be resummed to all orders in a s [19,20] to improve 
our prediction. 

Soft logarithms in the initial condition of the perturbative fragmentation function 
are process independent. We can hence resum them with NLL accuracy using the 
result presented in [6], which we do not report here for the sake of brevity. We 
present instead the results for the NLL resummation of process-dependent soft-gluon 
contributions in the MS coefficient function. 

In order to resum the large terms in eq. (3.1), we follow standard techniques 
[19], evaluate the amplitude of the process in eq. (2.1) at 0(as) in the eikonal 
approximation and exponentiate the result. The eikonal current reads: 



\J(Pt,Pb,Pg)f 



mi 



-2 



(Pt ■ Pb) 



(PfPg) 2 (Pt ■ Pg)(Pb ■ Pg) 



(3.3) 



For the sake of comparison with [19], we express the O(as) width in the soft ap- 
proximation as an integral over the variables 2 q 2 = (pt, + p g ) 2 x g and z = 1 — x g , 
with < q 2 < m 2 t {l - w) 2 {l - zf and < z < 1. The limits z -> 1 and q 2 -> 
correspond to soft and collinear emission respectively. In soft approximation, z ~ x b , 
the 6-quark energy fraction. We obtain: 

.jv-i 1 



T N (m t ,m w , Hf) 



C f 



TT JO 



dz- 



m?(i-™) 2 (i- 2 ) 2 dq 2 



-a s 



(1 - w) 2 (l - z) 2 Jo 



m t 2 (l- W ) 2 (l^) 2 



dq 2 as 



(3.4) 



In eq. (3.4) we have regularized the collinear singularity setting the cutoff q 2 > fx 2 F . At 
NLL accuracy level, this is equivalent to MS subtraction in dimensional regularization 
[6,21]. 

2 We point out that our definition of the integration variable q 2 is analogous to the quantity 
(1 — z)k 2 to which the authors of rcf. [19] set the scale for as for soft-gluon resummation in 
Drcll-Yan and Deep-Inelastic-Scattering processes. For small-angle radiation, q 2 ~ q 2 -,, the gluon 
transverse momentum with respect to the &-quark line. The variable z is analogous to z = 1 — E g /E q 
of ref. [19]. 
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In order to perform soft-gluon resummation to NLL accuracy a number of opera- 
tions have to be performed on this expression. We set the argument of as equal to q 2 
and, as far as the collinear-divergent term is concerned, we perform the replacement 



C F a s (q 2 ) A[a s (q 2 )} 



71 



(3.5) 



where the function A(as) was introduced in [19] and is detailed below. Moreover, 
the integral over q 2 of the non-collinear divergent term can be written, up to terms 
beyond NLL accuracy, as 



m 2 (l — w) 2 (l — z) 2 Jo 



m 2 (l- W ) 2 (l-*) 2 



dq 2 as(q 2 ) = as ( m t (1 — w ) 2 (l 



-z 



(3.6) 



This term describes soft radiation at large-angle, i.e. not collinear enhanced, and it 
is characteristic of processes where a heavy quark (the top quark in our case, the 
bottom quark in [15-18]) is present. It can be generalized to all orders by replacing 
eq. (3.6) according to: 



C 



__£ as ( m 2 (l _ w) 2 {1 _ z) 2^ S [ ft5 ( m 2 (l _ w) 2 {1 _ z) 2^ 



(3.7) 



This function is called T{a s ) in [15], B{a s ) in [16], S{a s ) in [17], D{a s ) in [18]. We 
now insert Eqs. (3.5-3.7) into eq. (3.4) and exponentiate the result. We obtain: 



In 



_ 1 [ f m1{l-wf{l-zf dq 2 



,1 Z N-1 _ 1 ( .m] 

A N = / dz— <^ / 

Jo i — z yjn 2 F 

+ S[a s (m 2 (l-w) 2 (l-z) 2 ) 



(3.8) 



We would like to evaluate eq. (3.8) to NLL level. The function A(as) can be expanded 
as follows: 



Mas) = E 



n=l 



7T 



(3.9) 



The first two coefficients are needed at NLL level and are given by [19,22]: 

= C F , (3.10) 



A^ = \c t 



V 18 6 , 
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(3.11) 



where Cf = 4/3, Ca = 3 and rif is the number of quark flavours, which we shall 
take equal to five for bottom production. 

The function S(as) can be expanded according to: 



S(a s ) = t(-) n S {n) 

n=l V 71 J 



(3.12) 
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At NLL level, we are just interested in the first term of the above expansion, which 
is given by: 

S<U = -C F . (3.13) 

The integral in eq. (3.8) can be performed, up to NLL accuracy, by making the 
following replacement [19]: 



- 1 -> -e i - 



-IE 



N 



— z 



(3.14) 



being the Heaviside step function. This leads to writing the following result for 
the function A^: 



A N (m t ,m w ,a s (n 2 ), fi, fi F ) = exp \nNg {1 \\) + g {2 \\ fx, li f ) 



with 



A = b a s {fi 2 ) In N , 
and the functions gW and g^ given by 



g {2) (A, //, n F ) 



+ 



2tt6 A 

AW 

2nb 

1 



[2A + (l-2A)ln(l-2A)] , 



(3.15) 
(3.16) 

(3.17) 



In 



m 2 (l — w) 2 

Ll 2 F 



-2 lE 



ln(l - 2A) 



4A + 21n(l-2A) + ln 2 (l-2A) 



+ 



2tt6 
27r6 n 



[2A + ln(l - 2A)] 
ln(l -2A). 



7T&o 



+ AW In 4- 



(3.18) 



In eq. (3.15) the term In NgW(X) accounts for the resummation of leading logarithms 
ag\n n+1 N in the Sudakov exponent, while the function g^ 2 \\, fi, fx F ) resums NLL 



terms a n s ln n N. 



Furthermore, we follow ref. [6] and in our final Sudakov-resummed coefficient 
function we also include the constant terms of eq. (3.2): 



t s N {m t) m W) a s {Li 2 ),Li,Li F ) = 



, . «s(/i 2 ) C F 

1 H K(m t , m w , fJ>F) 



2tt 



x exp 



\nNgW(\)+gW{\,»,» F j\. (3.19) 



One can check that the 0(cts) expansion of eq.(3.19) yields eq. (3.1). 

We now match the resummed coefficient function to the exact first-order result, 
so that also 1/JV suppressed terms, which are important in the region Xb < 1, are 
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taken into account. We adopt the same matching prescription as in [6]: we add 
the resummed result to the exact coefficient function and, in order to avoid double 
counting, we subtract what they have in common, i.e. the up-to-0(as) terms in the 
expansion of eq. (3.19). Our final result for the resummed coefficient function reads 3 : 

f^ s (m t ,m w ,a 5 (/i 2 ),/i,/i F ) = T s N {m u m w , a s (fi 2 ), fi, fi F ) 

- f^(mi,m H /,a 5 (/i 2 ),/i,/i F ) 

L J as 

T N {m t ,m w ,OLs{n 2 ),n,HF) , (3.20) 



+ 



where [f %\ as and [f w]a s are respectively the expansion of eq. (3.19) up to 0(as) and 
the full fixed-order top-decay coefficient function at 0(as), evaluated in Appendix 
B of ref. [14]. 

Before closing this section we would like to add more comments on the compar- 
ison of our resummed expression with other similar results obtained in heavy quark 
decay processes [15-18]. 

Besides the obvious replacement of a bottom quark with a top in the initial state, 
our work presents other essential differences. We have resummed large collinear 
logarithms as ln(m 2 /m 2 ), while Refs. [15-18] just address the decay of heavy quarks 
into massless quarks. Moreover, this work still differs in a critical issue. Those papers 
are concerned with observing the lepton produced by the W decay or the photon in 
the b — * X s 7 process, while we wish instead to observe the outgoing b quark. This 
is immediately clear from the choice of the x variable whose x — > 1 endpoint leads 
to the Sudakov logarithms. In our case it is the normalized energy fraction of the 
outgoing bottom quark; in [15-18] it is instead related to the energy of either the 
lepton or the radiated photon. 

The most evident effect of this different perspective is that an additional scale, 
namely the invariant mass of the recoiling hadronic jet, enters the results [15-18], but 
it is instead absent in our case. An additional function (called 7(0:5) in [15,16], C(as) 
in [17], B(as) in [18]) appears in those papers. The argument of as in this function 
is related to the invariant mass of the unobserved final state jet constituted by the 
outgoing quark and the gluon(s). It is worth noting that an identical function, called 
B [as (Q 2 (l — z))], also appears in the e + e~ [6] and DIS [19] massless coefficient 
functions, where it is again associated with the invariant mass of the unobserved jet. 
We do not have instead any B(as) contribution in our result. In fact, this function 
contains collinear radiation associated with an undetected quark, which we do not 
have since the b quark is observed. 

We also observe that to order a s the coefficient coincides with the corre- 
sponding of the function H [a s {ml(l — z) 2 )], which resums soft terms in the 



3 Alternative ways of matching, identical up to order as and differing in higher-order, subleading 
terms, are of course possible. 
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initial condition of the perturbative fragmentation function [6]. It will be very inter- 
esting to compare the functions S(as) and H(as) at higher orders as well. 

One final comment we wish to make is that, as expected, in our final result, 
eq.(2.5), which accounts for NLL soft resummation in both the coefficient function 
and the initial condition of the perturbative fragmentation function, ln n+1 iV terms 
do not appear, since they are due to soft and collinear radiation. Both the quarks 
being heavy, only the former leads to a logarithmic enhancement. Double logarithms 
are generated by a mismatch in the lower and upper q 2 integration limit over the 
Afas^g 2 )] function in the exponent of the resummation expression. In our case both 
of them have the same functional dependence with respect to z, i.e. (1 — z) 2 (see 
eq. (3.8) and eq. (69) of ref. [6]). The cancellation of the Q^ln"^ 1 ^ term can be 
explicitly seen at order as by comparing the large- N limit for the coefficient function, 
eq. (3.1), and the initial condition (eq. (45) of ref. [6]): the In 2 N terms have identical 
coefficients and opposite signs. 

4. Energy spectrum of the b quark 

In this section we present results for the 6-quark energy distribution in top decay. The 
6-quark spectrum in iV-space T N (m t , mw, I^f) is given by eq. (2.5). In the following 
we shall normalize IV to the full NLO width T [23], so that I\ = 1 will always hold. 
Results in avspace will be obtained by inverting numerically eq. (2.5) via contour 
integration in the complex plane, using the minimal prescription [24] to avoid the 
Landau pole. 

In order to estimate the effect of the NLL soft-gluon resummation, we compare 
our result with ref. [14] and use the same values for the parameters: m t = 175 GeV, 
m b = 5 GeV, m w = 80 GeV and A (5) = 200 MeV. 

In figure 1 we present the x b distribution according to the approach of perturba- 
tive fragmentation, with and without NLL soft-gluon resummation. For the scales 
appearing in Eqs. (2.6), (2.7) and (3.1) we have set hf = H = rn t and /io = Hof = n^b- 
We note that the two distributions agree for x b < 0.8, while for larger x b values the 
resummation of large terms x b — > 1 smoothens out the distribution, which exhibits 
the Sudakov peak. Both distributions become negative for x b — > and x b — > 1. As 
discussed in [14], the negative behaviour at small x b can be related to the presence 
of unresummed a^lna^ terms in the coefficient function. At large x b , we approach 
instead the non-perturbative region, and resumming leading and next-to-leading log- 
arithms is still not sufficient to correctly describe the spectrum for x b close to 1. In 
fact, the range of reliability of the perturbative calculation has been estimated to be 
x b < 1 - A/m b ~ 0.95 [6]. 

It is interesting to investigate the dependence of phenomenological distributions 
on the renormalization and factorization scales which enter the coefficient function (/x 
and hf) and the initial condition of the perturbative fragmentation function (/i and 
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Figure 1: fe-quark energy distribution in top decay according to the perturbative frag- 
mentation approach, with (solid line) and without (dashes) NLL soft-gluon resummation. 
In the inset figure, we show the same curves on a logarithmic scale, for x\, > 0.8. We have 
set jiF = A* = m t an d fJ-oF = fJ-o = mt- 



Hqf)- m particular, it is worth comparing the 6-energy spectra with and without soft 
resummation. For the scales /i and /xp we consider the values m t /2, mt and 2m<; for 
Hof and /x the choices are m&/2, m& and 2ra&. Figures 2 and 3 show the dependence 
of the Xb spectrum on the factorization scales /ip and /x f; the dependence on the 
renormalization scales \i and \i§ is exhibited in figures 4 and 5. 

We note that all distributions which include soft-gluon resummation exhibit a 
reduced dependence on the factorization and renormalization scales. 

Figure 2 shows that curves obtained using different values of /ip are almost 
indistinguishable once soft resummation is included; the unresummed plots exhibit 
instead a stronger effect of the chosen value for fip. A similar result also holds for 
the scale nop: the dependence of the plots on its actual value for Xb > 0.8 is small 
if soft logarithms are resummed and quite strong if the prediction is unresummed 
(figure 3). 

The choice of the value for the renormalization scale \x appearing in eq. (3.1) 
affects only the neighbourhood of the Sudakov peak of the resummed predictions, at 
^b- values very close to one (figure 4), where, as we have pointed out, our perturbative 
approach is anyway unreliable. The effect of the choice of the renormalization scale 
fi on the soft-resummed spectra is slightly larger than the one of \i and visible 
at x h < 1 as well (figure 5). As for the non-soft-resummed predictions, although 
all dashed curves in figures 4 and 5 seem to converge to same point for Xb — > 1, 
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0.0 0.2 0.4 0.6 0.8 1.0 



Figure 2: 6-quark energy spectrum for different values of the factorization scale /j,f, with 
(solid) and without (dashes) NLL soft-gluon resummation. The other scales are fixed at 
H = rrit, Ho = Hof = n^b- As in figure 1, in the inset figure, we present the same plots for 
large values of xi,, on a logarithmic scale. 



Solid: soft resummed 
Dashes: non resummed 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 3: As in figure 2, but for different values of hof- The other scales are fixed at 
fi = n F = m t , Ho = m b . 

the overall dependence on ji and jio for Xb < 1 is stronger than for the resummed 
predictions. 
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Figure 4: As in figure 2, but for different values of the renormalization scale \x. The other 
scales are fixed at \ip = m t , no = (J-of = n^b- 




Figure 5: As in figure 2, but for different values of the renormalization scale fiQ. The 
other scales are fixed at /x = /j,f = rnt, Hqf = mb- 



As a whole, one can say that the implementation of NLL soft-gluon resum- 
mation, along with the NLL DGLAP evolution for the perturbative fragmentation 
function, yields a remarkable improvement of our phenomenological results, since the 
reduced dependence on the choice of factorization and renormalization scales in the 
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region where the perturbative approach is reliable corresponds to a reduction of the 
theoretical uncertainty. 



5. Energy spectrum of 6-flavoured hadrons in top decay 

In this section we consider the problem of including a non-perturbative component 
on top of the perturbative result, so as to make predictions for observable 6-flavoured 
hadrons (like B mesons) in top decay. At the same time, we also account for the 
inclusion of NLL soft and collinear resummation. 

We write the normalized rate for the production of 6-hadrons B as a convolution 
of the rate for the production of b quarks in top decay, given by eq. (2.5), and a non 
perturbative fragmentation function D np (x): 

1 dT B , 1 r 1 dz dT b 



, K x B ,m t ,m w ,m b ) = ^ / — ^—(2, m t , m w , m b )D np (— \ (5.1) 
T dx B T Jx B z dz \ z J 

where x B is the B normalized energy fraction: 

xb = — , (5.2) 

1 — w mf 

Pb being the B four-momentum. Since D np (x) contains non-perturbative informa- 
tion, it cannot - for the time being - be calculated from first principles in QCD, but 
can only be extracted from data. We shall assume a universality property for such 
a function, and extract it from fits to ^-production data collected at LEP in e + e~ 
collisions. In particular, we can choose different functional forms for D np (x), and 
tune these hadronization models to the data. We shall consider a power law with 
two tunable parameters: 

g "^a,/ ?) = g(/; + 1 1 |a + 1) (l-,)V, (5.3) 

the model of Kartvelishvili et al. [30] 4 : 

D np {x; 5) = (1 + 5){2 + 8)(1- x)x 5 , (5.4) 
and the Peterson et al. model [31]: 

D np (x; e) = , A „ — . (5.5) 

v ' ; x[l - l/x - e/(l - x)] 2 K J 

In eq. (5.3), B(x, y) is the Euler Beta function; in (5.5) A is a normalization constant. 

4 We correct a typing mistake of ref. [14], where the normalization factor is the inverse of the cor- 
rect one. The numerical results of ref. [14] are nonetheless obtained using the correct normalization 
of the Kartvelishvili non-perturbative fragmentation function. 
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In order for our procedure to be self-consistent, care must be taken to employ the 
same underlying perturbative description in both the e + e~ — > bb process (where the 
non-perturbative contribution is fitted) and the t — > bW one (where it is used) . This 
will be ensured by using in both cases NLO, MS coefficient functions, along with 
a fully NLL soft-gluon resummed description, with the large collinear logarithms 
resummed to NLL accuracy by DGLAP evolution. For the coefficient functions in 
e + e~ annihilation we shall refer to [25]. 

Fits to data points can be performed either in x B -space, or, as recently advocated 
[10], in the conjugated moment space. When fitting in xb space we discard data 
points close to = and xb — 1 and consider ALEPH [26] data 5 in the range 
0.18 <xb< 0.94. The results of our fits are shown in table 1. 

The best-fit values for the parameters 
of the hadronization models are quite dif- 
ferent from the ones quoted in [14], where 
soft-gluon resummation in the e + e~ coeffi- 
cient function has not been used in the fits. 
The models in Eqs. (5.3) and (5.4) yield 
very good fits to the data, 6 as already found 
in [14], while the model (5.5) is marginally 
consistent. 

Using the results in table 1 we can give 
predictions for the spectrum of 6-flavoured 
hadrons in top decay. As in ref. [14], to ac- 
count for the errors on the best-fit param- 
eters, we shall plot bands which correspond 
to predictions at one-standard- deviation con 
fidence level. In figure 6 we show our pre- 
dictions for the Xb distribution using the 
three models fitted to the ALEPH data. 
At one-standard-deviation confidence level, 
the three predictions are different, with the 

Peterson model yielding a distribution which lies quite far from the other two and 
peaked at larger x^-values. Within two standard deviations, the predictions obtained 
using the models (5.3) and (5.4) are nonetheless in agreement 7 . We also note that 



a 


0.51 ±0.15 


(3 


13.35 ±1.46 




2.56/14 


5 


17.76 ±0.62 


X 2 (5)/dof 


10.54/15 


e 


(1.77 ±0.16) x 10~ 3 


X 2 (e)/dof 


29.83/15 



Table 1: Results of fits to e+e~ -> bb 
ALEPH data, using matched coefficient 
function and initial condition, with NLL 
DGLAP evolution and NLL soft-gluon re- 
summation. We set A( 5 ) = 200 MeV, 
Hof = A*o = m b = 5 GeV and fip = 
[i = \fs = 91.2 GeV. a and f3 are the pa- 
rameters in the power law (5.3), 5 refers 
to (5.4), e to (5.5). The fits have been 
performed neglecting the correlations be- 
tween the data points. 



5 Good-quality data are also available from SLD [27] on b- flavoured mesons and baryons. While 
this paper was being finalized, the OPAL [28] and DELPHI [29] Collaborations also published new 
results, which are fully compatible with the ones from ALEPH. 

6 We have also fitted the SLD data and found qualitatively similar results. However, the value 
for the parameters which best fit the SLD data are different from the ones obtained for ALEPH 
and quoted in table 1. 

7 Thc differences between the various models mainly originate from the varying quality of the 
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figure 6 is qualitatively similar to figure 4 of ref . [14] , where soft-gluon resummation 
in the coefficient function had not been implemented. In fact, the different pertur- 
bative content of these curves is now compensated by different values for the non- 
perturbative parameters a, j3, 8 and e, accordingly set by the fitting procedure. This 
hadron-level similarity does however not lessen the importance of the higher degree of 
reliability of the perturbative-level calculation provided for by the soft-gluon resum- 
mation: Such an accurate calculation can be used as a firmer starting point for testing 
and fitting non-perturbative models. On the other hand, the inclusion of perturba- 
tive resummation cannot be expected to improve the agreement between hadron-level 
results obtained with different phenomenological non-perturbative ansetze. 




Figure 6: xb spectrum in top decay, with the hadronization modeled according to a power 
law (solid lines), the Kartvelishvili et al. (dashes) and the Peterson (dots) model, with the 
relevant parameters fitted to the ALEPH data. The plotted curves are the edges of bands 
at one-standard-deviation confidence level. NLL soft-gluon resummation is included. We 
set [j>f = fi = rcii and [i$f = Mo = m b- 

An alternative, and probably better, way of determining and including non- 
perturbative information makes use of moment space perturbative predictions and 
data [10]. The full hadron-level result can be written in A^-space as the product of a 
perturbative and a non-perturbative contribution, = T b N D^f . For each value of A" 
one can then extract the corresponding D^f value from e + e~ data, with no reference 
whatsoever to a specific hadronization model, and use it to predict the same moment 
in top decay. The DELPHI Collaboration [29] has recently published preliminary 
results for the moments of _B-meson fragmentation in e + e~ collisions up to A^ = 5. 
From these data, and using the moments of the e + e~ perturbative contributions [6], 

fits to e+e~ data where, as the x 2 values in table 1 seem to suggest, a given model is sometimes 
not really able to describe the data properly, due to its too restrictive functional form. 
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one can extract D^f. The corresponding values can then be calculated making 
use of the results for T b N obtained in this paper. Calling a§ and a b N the moments for 
the production rate of B mesons (measured) and b quarks (calculated in perturbative 
QCD) in e + e~ annihilation, we have a B = cr^D^f and hence 

rr B 

N N N — 1 N b ' v°' D / 

Table 2 shows a practical implementation of this procedure. Predictions for the 
moments of _B-meson spectra in top decay are given, making use of the DELPHI 
experimental data. Two sets of perturbative results ([A] and [B]) are shown, the first 
using A^ 5 ) = 0.226 GeV and = 4.75 GeV, the second using the default parameters 
of this paper. As expected, the perturbative calculations and the corresponding non- 
perturbative components differ, but the final predictions for the physical results T^- 
are to a large extent identical. 





(x) 


<* 2 ) 


<z 3 > 


<z 4 > 


e + e data a B 


0.7153±0.0052 


0.5401±0.0064 


0.4236±0.0065 


0.3406±0.0064 


e+e" NLL a b N [A] 


0.7666 


0.6239 


0.5246 


0.4502 


e+e" NLL a b N [B] 


0.7801 


0.6436 


0.5479 


0.4755 


D% [A] 


0.9331 


0.8657 


0.8075 


0.7566 


D% [B] 


0.9169 


0.8392 


0.7731 


0.7163 


t-decay NLL T b N [A] 


0.7750 


0.6417 


0.5498 


0.4807 


t-decay NLL T b N [B] 


0.7884 


0.6617 


0.5737 


0.5072 


t-decay [A] 


0.7231 


0.5555 


0.4440 


0.3637 


t-decay L^ [B] 


0.7228 


0.5553 


0.4435 


0.3633 



Table 2: Experimental data for the moments a B from DELPHI [29], the resummed e + e~ 
perturbative calculations for a b N [6], the extracted non-perturbative contribution D^f. Us- 
ing the perturbative results T b N , a prediction for the physical observable moments L^ is 
given. Set [A]: A( 5 ) = 0.226 GeV and m b = 4.75 GeV, set [B]: A^ 5 ) = 0.2 GeV and 
= 5 GeV. The experimental error should of course be propagated to the final predic- 
tion. 



6. Conclusions 

We have discussed soft-gluon resummation in top quark decay within the framework 
of perturbative fragmentation which, by making use of the DGLAP evolution equa- 
tions, allows one to resum with NLL accuracy the large logarithms as In (mf/m 2 ) 
which appear in the 0{ag) massive calculation. The MS coefficient function and 
the initial condition of the perturbative fragmentation function contain terms which 
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become arbitrarily large for soft-gluon radiation. Such soft terms are process depen- 
dent in the coefficient function and process independent in the initial condition of 
the perturbative fragmentation function. 

We have performed the resummation of soft-gluon effects in the coefficient func- 
tion of the top-to-massless-quark decay process with NLL accuracy, and matched the 
resummed result to the full 0(a s ) one. This result has then been combined with the 
one of ref . [6] , which resums NLL process-independent soft-gluon contributions in the 
initial condition of the perturbative fragmentation function, to produce a resummed 
prediction for the full top-to-massive-bottom decay. 

We have presented the resummed 6-quark energy spectrum in top-quark decay 
and compared it with the unresummed prediction. We have found that at large Xb 
the implementation of soft-gluon resummation has a visible impact: The Xb spectrum 
is smoothed out and shows the characteristic Sudakov peak. Our prediction for the 
6-energy spectrum is negative for — > and x b — > 1; we have interpreted the 
behaviour at small Xb as due to unresummed liaxb contributions in the coefficient 
function and at large Xb to missing non-perturbative effects. 

We have investigated how the 6-energy spectrum varies if we choose different 
values for the factorization and renormalization scales which enter our calculation. 
We have found that after including NLL soft-gluon resummation the distributions 
exhibit very little dependence on the choice of such scales, which corresponds to a 
reduction of the theoretical uncertainty of our prediction. 

Finally, we have made use of ^-flavoured hadron data in e + e~ collisions to extract 
a non-perturbative contribution, and we used it to calculate hadron-level predictions 
in the top decay process. This has been done both in space, tuning various 
hadronization models to e + e~ distributions, and in N space, directly extracting the 
values of a few moments of the non-perturbative fragmentation function. In both 
cases care has been taken to make the procedure self-consistent by employing the 
same perturbative description in both the e + e~ and top decay processes. 
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